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Computing the perturbation spectrum in the recently proposed Island Cosmology remains an open 
problem. In this paper we present a classical computation of the perturbations generated in this 
scenario by assuming that the NEC-violating field behaves as a classical phantom field. Using an 
exactly-solvable potential, we show that the model generates a scale-invariant spectrum of scalar 
perturbations, as well as a scale-invariant spectrum of gravitational waves. The scalar perturbations 
can have sufficient amplitude to seed cosmological structure, while the gravitational waves have a 
vastly diminished amplitude. 



I. INTRODUCTION 

Recent cosmological obervations [lj seem to be con- 
sistent with the hypothesis that we live in a flat A- 
dominated Universe seeded by primordial fluctuations 
that were predominantly adiabatic and nearly scale- 
invariant. The most successful class of models leading 
to such a Universe are the Inflationary models (see, 
for instance 0, 0, and @ for recent reviews). 

Recently a new cosmological model has been proposed 
in Ref. Q called Island Cosmology. In this scenario, the 
presently observed state of the Universe (i.e., inflating 
with a cosmological constant), is considered to be the 
eternal state. Large-scale explosive events which involve 
local violations of the Null Energy Condition (NEC) lead 
to the formation of "islands" of matter, one of which 
is our observed universe. The cause of these explosive 
events is attributed to quantum fluctuations in some 
field. While a complete calculation of density pertur- 
bations in this scenario (which would require a proper 
treatment of the back reaction in quantum gravity) is not 
attempted, the authors show that any other free scalar 
field responding to the metric fluctuation would produce 
a scale-invariant spectrum. 

In this paper we consider a special version of Island 
Cosmology where we model the behavior of matter dur- 
ing the NEC-violating event as a classical phantom field 
in order to compute the perturbation spectrum. The dif- 
ferent stages of our model are described below (see Fig. 
®: 

1. The de Sitter Phase: We start with the assumption 
that the initial state of the Universe is de Sitter 
space inflating with the currently observed value of 
the Hubble constant. 

2. The Phantom Phase: A quantum fluctuation in 
some field (chosen to be a light scalar field for sim- 
plicity) in the expanding phase of de-Sitter space- 
time, and occuring over a Horizon-sized volume 
drives the Hubble constant to a high value. Such 
a fluctuation necessarily violates the Null Energy 
Condition within this volume. In this treatment 
we make the working hypothesis that during this 
phase the energy content within a (current) Hubble 
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FIG. 1: Sketch of the behavior of the Hubble length scale 
with conformal time, n, in the Island model, and the evolu- 
tion of fluctuation modes. At early times, inflation is driven 
by the presently observed dark energy, assumed to be a cosmo- 
logical constant. As the cosmological constant is very small, 
the Hubble length scale is very large - of order the present 
horizon size. Exponential inflation in some horizon volume 
ends not due to the decay of the vacuum energy as in in- 
flationary scenarios but due to a quantum fluctuation in the 
time interval (rii,rjf) that violates the null energy condition 
(NEC). The NEC violating quantum fluctuation causes the 
Hubble length scale to decrease. After the fluctuation is over, 
the universe enters radiation dominated FRW expansion, and 
the Hubble length scale grows with time. The physical wave- 
length of a quantum fluctuation mode starts out less than 
H^ 1 at some early time r\i. The mode exits the cosmological 
horizon during the NEC violating fluctuation (?7exit) and then 
re-enters the horizon at some later epoch (gentry) during the 
FRW epoch (The modes are drawn as straight lines for illus- 
trative purposes only, they actually grow in proportion to the 
scale factor). 

sized volume of the Universe behaves like phantom 
energy, i.e., a classical perfect fluid with equation 
of state w = p/p < — 1, where p and p stand for 
pressure and density respectively, and p > 0. The 
Hubble length decreases sharply during this period, 
though the Universe keeps expanding. 

3. The FRW Phase: After the explosive event, the 
Hubble constant is large, and classical radiation fills 
the volume. Rapid interactions thermalize this ra- 
diation, and this part of the Universe then follows a 
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radiation-dominated FRW evolution leading to the 
Universe that is observed today. 

4. The Dilution Phase: Eventually cosmic expansion 
dilutes out this aggregation of matter, and the vol- 
ume of space is restored to its initial de Sitter state. 

In the remainder of this paper we will show that this 
model predicts a (nearly) scale invariant spectrum of 
scalar perturbations with amplitude sufficient to seed 
structure, as well as a scale invariant spectrum of ten- 
sor perturbations of considerably smaller amplitude. 

A similar situation was studied by Y.S. Piao in Ref. 
@ using a different calculational approach. 

The layout of the paper is as follows: Section [H] de- 
scribes the model in greater detail. Section lull discusses 
the assumptions and outlines our calculational strategy. 
Section IIVI presents a detailed calculation of both the 
scalar and tensor power spectra, and Section |V] summa- 
rizes our findings. 



II. THE MODEL 

In the following subsections, we describe in more detail 
the three stages of the model, stating the assumptions 
involved in each stage. 

We set our notation as follows: let us choose our cos- 
mic time coordinate t such that the phantom phase (de- 
scribed in Section lll Bl) begins at t = ti = 0. In this paper 
we find it more convenient to work in conformal time (rj) 
where dt — a(rj) drj and r\ ranges between (-co, 0) as t 
goes from — oo to +oo. We set our conformal time coor- 
dinates such that the period of phantom cosmology lasts 
between r\ — r/i = — 1/H& and r\ = r/f. The above choices 
of U and r]i are arbitrary and for convenience. Primes 
denote derivatives with respect to conformal time, and 
dots denote derivatives with respect to cosmic time. We 
adopt the convention that the suffix i denotes the value 
of a quantity at r\ = rji and the suffix / denotes the value 
at 1} = r)f. 



A. The de Sitter phase 

This phase represents the initial state of the Universe, 
before the onset of the phantom behavior. In this phase, 
we assume that the Universe is de Sitter space inflating 
due to the observed dark energy, which we assume is a 
cosmological constant. The Hubble parameter (H) has 
the same value that it has today, which we call H\ . This 
expanding de Sitter background can be part of a classical 
de Sitter spacetime with no beginning and no end, with 
early contraction and then expansion. We will only con- 
sider the expanding phase of the de Sitter spacetime in 
the following discussion. 

We assume that the matter content of the Universe is 
a classical scalar field (4>) having a Lagrangian C given 



by: 

and a stress tensor given by: 



Tpv = Xdfj, (4>)d u (4>) - g^C 



(1) 



(2) 



gt^f represent the components of the metric, and V((p) 
represents the potential. For the sake of generality, we 
have inserted the constant A which determines the sign 
of the kinetic term. Obviously, A = +1 for an ordinary 
classical scalar field, and A = —1 for a phantom field. 

The equation of motion of the field is the Klien Gordon 
equation: 

-g^dM + ^ = (3) 



The scale factor a(t) and Hubble value H(t) during 
this period can be written as follows: 
For (-00 <t<ti=0): 



a(t) = e H ^ 
Hit) = ^=H A 

In terms of conformal time, for (— oo < ij < rji): 

1 



a(rj) = - 
H(t]) 



r)H A 



a{rj) rj 



(4) 
(5) 



(6) 
(7) 



B. The Phantom Phase 

In this phase, lasting between the times rji and rjf, the 
Universe undergoes an NEC-violating quantum fluctua- 
tion. We model this phase by assuming that the matter 
content of the Universe behaves like a phantom field <p p 
for the duration of the fluctuation. 

We assume that this hypothetical phantom field <p p is 
classical, i.e., its Lagrangian and stress tensor are given 
by Eq. (J) and Eq. (J2J), and it satisfies the Klien Gordon 
equation Eq. ©. 

Using the above equations, one can readily determine 
the pressure p and energy density p of the field. These 
work out to: 



p = \Y,Tu =\^~a\ V )V{4> P ) (8) 
<b' 2 

P= T 00 =\^ + a 2 {f 1 )V{4> p ) (9) 

Clearly, p + p = X()> p indicating, as one would expect, 
that the NEC is violated during this period if our matter 
field is phantom (A = —1). 
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Also during this phase, the Hubble horizon size drops 
from H^ 1 to HJ 1 . We assume that this drop is linear 
in cosmic time, ending at time tf. The validity of this 
assumption, as well as its implications on the matter con- 
tent of the Universe are discussed later. Thus, 



The scale factor in this radiation-dominated epoch can 
therefore be written as 



a(t) = a f , 



(18) 



H- X {t) = HJ 1 - at (for 0<t<tj 



(10) 



Here a is a dimensionless parameter measuring the rate 
at which the horizon size changes during the phantom 
phase. We assume that the quantum fluctuation is very 
abrupt, and hence a is very large. 

Using the definition of conformal time, it is easy to 
deduce that during this phase (rji < r) < t] f) , 



a(rj) = 
H( V ) = 

H(rj) = 



[-a-H A (l 

<m 

a(v) [-a 
H{rj) 



a)rf\ i+° 



a{rj) 



[-a - H A (1 + a)n] 
Ha HvT 



(11) 
(12) 

(13) 



We also need to address the nature of the back- reaction 
of matter on geometry. Let us make the working hy- 
pothesis that the back-reaction is fully described by the 
Friedmann equation. 

The actual time duration of this phase can be calcu- 
lated by demanding continuity of the Hubble value at 
t = tf (or equivalently, 7/ = r/f). Thus we require (see 
Eq. 1131)), 



H(Vf) = H f 
Solving the above equation for 77/, we obtain 

1+-" 

' / " A % > 

Vf 



H A (l + a) 



(Ha 



(14) 



(15) 



We assume that Hf » Ha and since a — > 00, the first 
term in the square brackets can be ignored leaving us 



Vf 



H A (1 



O 



1 



l + a)Hf 



(16) 



From this we can compute the duration of the phantom 
phase in conformal time (Ary), as follows: 



(17) 



Aif = rif - m 

= -(— 

H A \a+l 

Again, since a — > 00, this is a vanishingly small interval. 

C. The Radiation Dominated FRW Phase 

In this epoch we have a volume of space of Hubble 
length HJ 1 filled with classical radiation. Rapid inter- 
actions thermalize the radiation, after which this volume 
follows a standard FRW evolution. 



In terms of the conformal time, this reduces to 

,2 



a(rj) = a f +a f H f (r)-r)f) 



III. CALCULATIONAL STRATEGY 



(19) 



We make two key assumptions to facilitate our calcu- 
lation, which we discuss below. These are: 

1. During the NEC- violating explosive event the en- 
ergy content of the Universe behaves as a phantom 
field. 

2. During the NEC violation, the drop in the Hubble 
scale is linear in cosmic time. 

Assumption Q is an attempt to model the behavior 
of the matter field during the NEC-violating event. To 
calculate density fluctuations due to fluctuations in the 
NEC-violating field, one needs a suitable model for the 
evolution of the field itself during the NEC- violating fluc- 
tuation. This evolution is quantum and not described 
as a solution to some classical equation of motion. For 
the purpose of this calculation, we have made the sim- 
plifying assumption that the matter field behaves in the 
same manner as a classical object that would also violate 
the NEC and produce the same effect on the spacetime. 
Of course this purely classical treatment cannot substi- 
tute for a rigorous quantum mechanical treatment of the 
NEC-violation, but we hope that it captures the essential 
elements of the physics involved. 

Assumption (J2J) can be justified considering that the 
drop in the Hubble length need not be linear throughout 
the explosive event, but only during the window of time 
5t that it takes for the scales observed today to leave the 
horizon. Since the fluctuation itself is very short lived, 
the drop in H~ l (t) over St can be well approximated to 
be linear. 

We now turn our attention to computing the spectrum 
of perturbations that would be generated in this cosmo- 
logical model. Our plan of action is the following: 

1. Working in k (momentum or wavenumber) space, 
we first find expressions for the Mukhanov variable 
Vk(rf) in all the three stages of the model. The 
Mukhanov variable is a gauge- invariant linear com- 
bination of matter and metric perturbations repre- 
senting the true dynamical degrees of freedom of 
the system. It is fully described in 0. 

2. The unknown coefficients that arise in the above ex- 
pressions are then determined by demanding conti- 
nuity of Vk (j]) and its time derivative v' k (77) at tran- 
sition times r\i and iff. 
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3. The adiabatic density perturbation responsible for 
structure in the Universe is conveniently character- 
ized by the curvature perturbation 1Z seen by co- 
moving observers. Once we fully determine Vk(r]) 
in the radiation dominated phase, we obtain the co- 
moving curvature perturbation spectrum from the 
relation (see, for example, 0) 



k 3 



z ( r l) FRW 

z being defined by the relation 



H 



IV. THE PERTURBATION SPECTRUM 



In the next three subsections l|IV Al llV Bl and lTV ( "(I . we 

find expressions for in the three stages of the model. 
The unknown coefficients are determined in IIVDI and 
the scalar and tensor power spectra are computed in llV El 
and IIV Fl respectively. 



A. Vk in the de Sitter Phase 

For a scalar field in de Sitter space, the gauge- invariant 
Mukhanov variable Vkirj) can be shown to satisfy the 
equation (see, for instance, 



a"(v) 
a(rj) 



(22) 



The solutions to this equation are the Bunch-Davies 
mode functions: 



-ihrj 



I 
kl] 



(23) 



B. Vk in the Phantom Phase 



Given our choice of gauge, the metric perturbation 
$(77,2?) coincides with the gauge-invariant Bardeen po- 
tential (see, for example 9]). 

The phantom matter field (j} p (r],x) is perturbed as fol- 
lows: 



(f> p (r),£) = <t>(rj) + 5(j)(x,rj) 



(25) 



(20) 

2. Evolution of the perturbations 

To find time evolution of the perturbations, we use the 
perturbed Einstein equations up to first order. The i-i 
(21) component of the zero-th order equations reads: 



2— = 8ttG 
a 



A^-a 2 U(0) 



while the 0-0 component reads: 



Aa 2 + 3 ( — ) = 8vrG 



f/2 



(26) 



(27) 



Adding these equations, one obtains the familiar rela- 
tionship: 



U 2 -W = 4nG\<f> K 



(28) 



The i-i, 0-0 and 0-i components of the first order equa- 
tions are respectively: 



Aa 2 $ - 2$ 



4$. 



a 



8ttG 



a 2 V{<f>)<P - ^a 2 U (50 + ^\(/)'5<f>' - A$0' : 



Aa 2 $-3— $' + V 2 $ = 
a 



8ttG 



a 2 U(</>)$ + -a%<50+-A<^ 



$' + —$ = ±1tG\<j>'t 



(29) 

(30) 
(31) 



In this case we will calculate starting from first prin- 
ciples. 



1. Matter and metric perturbations 

The first step is to perturb the matter and met- 
ric. Working in longitudinal gauge and assuming no 
anisotropic stress, the scalar metric perturbations are 
written as: 

1 + 2$ 

2 i -1 + 2$ 
~ " I -1 + 2$ 

-1 + 2$ 



(24) 



(where = dVj d(f>) . 

Perturbing the Klein-Gordon equation Eq. © using 
Eq. 125|) yields, at zero-th order: 



A [</>" + 2H(t>'} = -a 2 V<p 



(32) 



Using Eq. (ggjl. Eq. Pty. Eq. ffOjl. Eq. and Eq. 

one obtains the equation of motion of $(ry,a;): 



$"-V 2 $+2 [H - $'+2 H' - H-f- $ = (33) 



Applying the Fourier transform: 

d 3 k 



(2tt) 3 / 2 



$fc(?7)e 



ik.. 



(34) 
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we obtain: 



where xiv) satisfies the differential equation 



2 n 



2 In' -H- 



(35) 

Note that Eq. Ij35|l is independent of A, indicating that 
it has the same form for a phantom field as it would 
for a normal field. This is a surprising result, since all 
the equations used to derive Eq. (|35[1 are A dependent. 
Physically, this result implies that the evolution of the 
metric perturbation is insensitive to whether the matter 
content of the Universe is normal or phantom. 

To solve Eq. Ij35|l . we need to determine the dynamics 
of the phantom field (^(rj), which is in turn determined 
by the potential V(4>). We choose a particular form of 
the potential which allows for a solution in closed form: 



2 (3 + a) 



2aH A 
K 



(36) 



Here K is a constant that has the dimensions of mass 
squared, and sets the scale of the potential. 

For this potential, it is easy to verify that the exact 
form of <p(rj) which satisfies Eq. I|32ll (with of course A = 
— 1 as is the case in the phantom phase) is: 



K 



H A {l + a] 



In [-a - H A (1 + a)ri] (37) 



Further, to facilitate the back-reaction as discussed in 
Section 111 Til and satisfy our ansatz given by Eq. I|1(J|) 
we must require that our field satisfies the Friedmann 
equation Eq. (|2*7|) . The result of this is to fix the value 
of K: 



K = H 



3a 
4^G 



(38) 



The matter field <p{rj) has the interesting property that 



= (l + a)H( V ) 



(39) 



Using Eq. l|3l?)). Eq. (JSSJ reduces to 

$fc - 2Ha& k + [fc 2 - 2H 2 a] $ fe = (40) 

This is a familiar second order differential equation of the 
form: 



< + P(r ? K + Q(r / )$ fc = 



(41) 



with 



Q(v) 



-2H(r])a 

k 2 - 2H 2 {r])u 



which has the solution (see e.g. |Tl| 1 

$> k (v)=e-^ P W dr >x(v)=a a X (ri) 



(42) 



x"(v) 

In our case, this reduces to 

x"(v) 



\p'--P 2 
2 4 



k 



a + H A (l + a)i]} 



X(V) = (43) 



X(v) = (44) 



At this point it is convenient to temporarily switch to a 
new time variable x defined by 



x = a + H A (l + a)r) 



(45) 



Note that x = — 1 when r\ — r\i and from Eq. illtj|) . x = 
O [H\/Hf] or x ~ when r\ — r)f. 

With x as the time variable, Eq. (|44|l takes the form 



d 2 x{x) 
dx 2 



m — l p 



1\ 1 



4 / x 2 



where m and p are defined by 
k 



H A (l + a) 



1 



P - 



4 {a + lf 



1/2 



X (x) - (46) 

(47) 
(48) 



The solution to Eq. g^l can be written in terms of Bessel 
functions as 

x(x) = \Jmx [A m Jp(mx) + B m Y p (mx)] (49) 

where A m , B m are constants and J p and Y p denote Bessel 
functions of the first and second kind of order p respec- 
tively. 

Putting together Eq. and Eq. ljl2*)l . we conclude 
that the time evolution of is fully described by the 
equation 

®k(v) = a a ^/mx[A m J p (mx) + B m Y p (mx)} (50) 

where m, x and p are defined by Eq. (|47|l . Eq. 145J1 and 
Eq. l(4*%|l respectively, and a is defined in Eq. lfTT)|l . 



3. Calculating the Mukhanov Variable 

The Mukhanov variable Vk can also be defined by the 
relationship (see, for example 0): 



(51) 



Where IZk is the curvature perturbation on co-moving 
hyper-surfaces at length scale k. 

We are first going to compute IZk and then z and then 
use Eq. (|51|) to compute Vk- 
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a. Calculation of 7Z k ■' The curvature perturbation 
on co-moving hyper-surfaces (IZk) is defined by the rela- 
tion (in momentum space) Q: 



Noticing that for this space, we have 

W 1 



(52) 



(53) 



and using equations Eq. (|31[l and Eq. (|28|l . we can elim- 
inate the S<p in Eq. l|52[l . giving us an expression for TZ k 
involving and $' fc as the only first order variables: 



i 



(54) 



(Again, the absence of A indicates that the expression for 
lZ k is insensitive to whether the field is real or phantom.) 

b. Calculation of z: In the definition of z Eq. 121|) 
we eliminate </>' using Eq. (|28|) . and introduce a using 
Eq. to get 



(55) 




amp\ 



where mpi is the reduced Planck mass defined by mp\ — 
a/1/87tG, where G is Newton's gravitational constant. 

c. Final expression for v k : Combining Eq. (|51|l . 
Eq. f54^l and Eq. (|55|) . we can express Wfe entirely in terms 
of $fc and $1 as follows: 



v k = a 



-AttGX 



an 



(56) 



D. Calculation of Unknown Constants 

The above calculations produced four unknown con- 
stants A m , B m , in Eq. l|4^|) and a k and f3 k , in Eq. JSSJl. 
These constants can be determined by demanding conti- 
nuity of Vk and its time derivative at the two transition 
times r\i and 77/ . 

To determine A m and B m , we perform the above 
matching process at r\ — rji = — or in terms of x, 
(from Eq. (03) j &t x = Xi = —I . In other words, we 
need to simultaneously solve the equations 

v k (de Sitter) I ( ?7= _ 1 / Ha ) = v k (phantom) |( x= _i) 
v' k (de Sitter) I (^ = _i/jj a ) = ^(ph antom )l(a;=-i) (59) 

to find the unknowns A m and B m . 

The expression obtained for v k in the phantom phase 
(by substituting the values of the above constants) is 
fairly complicated. However, since we are only inter- 
ested in the super-horizon modes, we can make the ap- 
proximation that k (or m) — > 0, and use the appropri- 
ate asymptotic forms of the Bessel functions J p {mx) and 
Y p (mx). Also, since a is large, from Eq. p ~ ~. 

With these simplifications the expression for v k (in the 
phantom phase) reduces to: 



v k (x) 



v 'k( x ) 



ie lk / HA a ( - 1+a '>H A (-l + xa) 



\/2fc 3 /2(l + a ) 



+0 



ie ik/H Aa (l+ a ) H 2 

V2k 3 / 2 x + 



(60) 
(61) 



Now having fully determined the form of v k (77) during 
the phantom phase, we can determine the coefficients a k 
and /3 k in Eq. I|58|l to determine v k (rj) in the final FRW 
phase. In particular, we need to solve simultaneously the 
equations: 



C. v k in the FRW Phase 

Here we have a gauge field in a geometry described by 
Eq. H19|) . From the theory of cosmological perturbations, 
we know that v k satisfies 6] 



z"(r,) 



z{ri) 



v k (v) = 



(57) 



But since z(rf) = 2mpia(r]) in this space, the above equa- 
tion simply reduces to Eq. . Solving Eq. (|22l for this 
space, we find that v k is described by the solution: 



Vk(r) 



ikr , n —ikr 

a k e + (i k e 



(58) 



where a k and (3 k are constants of integration and r 
77 — i]f > 



u fe (phantom) | (a:::i0) = v k (FKW)\ {v=rif) 
v' k (phantom) | (a ,~ 0) = v' k (FRW)\ (v=flf) 



(62) 



For brevity, let us call the leading term in k in the ex- 
pression for v k (Eq. (fol^ at 77 = 77/ (or x ~ 0) as h and 
the leading term in k in the expression for v' k (Eq. (|61l) ') 
at x ~ as h- Thus we have 



and, 



h = 



^3/2(1 + a) 



lei k/H Aa (l+ a ) H 2 A 



(63) 



V2fc 3 /' 



x f 



ie i/kH Aa (.2+ 2a ) H 2 



V2fc 3 / 2 



(64) 
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Where in the last manipulation we have used the result 
that —Xf = a j ( 1+a \ which follows from the form of the 
scale factor during the phantom phase (Eq. Ulljl l and the 
definition of x (Eq. I]45[ll. Eq. l(6"21) now implies that 

a k =\[h + |] 

Using Eq. at (r) = r]f) to find 

aJH A = H f 

we note from equations Eq. (|64|) that, 

h _ 1 1 

h ~ (l + a)a f +a H A ~ a f H f (l + a) 

since both a and Hf are large. Eq. I|65|) now reduces to 

1 \h~ 

2 ik 
1 ' 

~2 



< 1 



(65) 
(66) 

(67) 



OLk 



0k 



(68) 



Hence the form of Vk{r]) in the phantom phase becomes: 



2ik 



2ik 



= h 



sin(ki 
kr 



(69) 



E. Determination of the Scalar Power Spectrum 

Now we are in a position to determine the power spec- 
trum of the co-moving curvature perturbation in the 
FRW space using Eq. I|20[l which gives 



k 3 
2it z 



sin(ki 



kr zfrw(?7) 



(70) 



ajHfT from 



Making the approximation apRw( r ?) 
Eq. (|19l) since at the time a mode re-enters the horizon 
■q ^> rif, using the relation Eq. (|S5|> . and taking the limit 
k — ► 0, we finally obtain 



Vn 



2^2 



ik/H A 



2V2k 3 / 2 m P i 



Hf 



\47rmpi J 



(71) 
(72) 



Hence we find that our model produces a (nearly) scale 
invariant spectrum of cosmological perturbations, with 
amplitude set by Hf/mp\. If we assume Hf ~ 10 14 GeV 
(approximately GUT scale), then the power spectrum 
matches the COBE DMR observations [ljj of CMB tem- 
perature fluctuations of order 10 -5 . In other words, the 
perturbation spectrum can have an amplitude sufficiently 
large to seed the cosmological structure that we see to- 
day. 



F. Determination of the Tensor Power Spectrum 

We know from the theory of cosmological perturba- 
tions (see, for example, |6j) that gravitational waves 
are essentially equivalent to two massless scalar fields 
(for each polarization) up to a renormalization factor of 
2/mp\. Hence we can write the tensor power spectrum 
Vtzt as 



Vtzt = 2 — —V-Rip 



(73) 



pi 



where the first factor on the right comes from the two 
polarization states, the second represents the renormal- 
ization mentioned above, and V-ji^ is the spectrum of 
perturbations of a massless scalar field ip other than, and 
not interacting with, the NEC- violating field. The V'-field 
perturbations can be computed using essentially the same 
machinery as above, with a small difference in the final 
step: 

1. We solve Eq. I|35(l with ip representing the matter 
field. 

2. We find expressions for Vk in the three stages of the 
model (up to constants of integration). 

3. We demand the continuity of Vk and v' k at rji and rjt 
to fully specify the expression for Vk in the FRW re- 
gion (that is, evaluate the undetermined constants 
obtained in the previous step), and from this, com- 
pute the perturbation spectrum Vn.^ using the re- 
lation 



2^2 



a(rj) 



(74) 



FRW 



The final expression for the power spectrum of tensor 
perturbations turns out to be 



H A 



( H A 



a 3 innp\ 



\7rmp1 



(75) 



The last step follows because it is easy to show (by Taylor 
expanding a{rf) about the point r\ = rji) that the scale 
factor hardly changes from its initial value of 1 during 
the NEC- violating event. 

This result agrees (up to 0(1) numerical factors) with 
the corresponding result obtained in pj . 

A similar scenario was investigated by Y.S Piao 
and while he obtains a scale invariant spectrum of scalar 
perturbations, the tensor perturbation spectrum in his 
calculation turns out to be blue-shifted. The source of the 
discrepancy could lie in his assumption that the canon- 
ical relationship satisfied by Vk in the case of the tensor 
spectrum has a time-dependent mass given by a"(r])/a(r]) 
(Eq. (12) in Q) during the phantom phase. In our ap- 
proach we derive Vk from first principles and find that 
the time dependent mass can have a more complicated 
form. 
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V. CONCLUSIONS 

To summarize, in this paper we have investigated a cos- 
mological model in which islands of matter are created 
from NEC-violating explosive events in a cosmological 
constant dominated inflating Universe. Our approach 
was to model the behavior of the matter field during the 
NEC-violating fluctuation as a classical phantom field. 
While an ideal approach would be a full quantum me- 
chanical calculation, it is hoped that our classical calcu- 
lation can capture the essential elements of the physics 
involved. 

Our calculations yield, (using an exactly solvable po- 
tential), an adiabatic spectrum of scale-invariant pertur- 
bations, whose amplitude is determined by the value of 
the Hubble constant at the end of the NEC- violating fluc- 
tuation. If we assume the latter to be approximately 
GUT scale, the perturbation spectrum turns out to have 
an amplitude sufficient to seed the cosmological structure 
seen today. We also obtain a scale-invariant spectrum of 



gravitational waves of amplitude set by (H\/mp\). 

The essential elements of Island Cosmology, namely, 
the A-sea and the scale invariant perturbation spec- 
tra, are supported by current observations. The NEC- 
violating quantum fluctuations which create the matter 
islands are a clear consequence of the basic formalism of 
quantum field theory. Island cosmology therefore holds 
out the tantalizing notion that the observed Universe 
could have had a purely quantum field theoretic origin. 
This paper lends further support to this idea. 
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